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Stress  Intensity  Solutions  for  Corner  and  Through-the-Thickness 
Cracks  at  a  U-Shaped  Notch 

JAMES  A.  HARTER 

Sustainable  Structures  Branch,  United  States  Air  Force  Research  Laboratory,  Wright-Patterson  AFB,  OH 


ABSTRACT 

Stress  intensity  factor  solutions  for  comer  and  through  cracks  growing  from  a  u-shaped  edge  notch  are 
required  to  improve  the  accuracy  of  damage  tolerance  analyses  for  continuing  damage  scenarios.  Closed- 
form  K-solutions  for  through-the-thickness  and  comer  cracks  under  remote  axial  loading  were  developed 
using  3-D  finite  element  (FE)  models  for  wide  plates.  Finite  width  corrections  were  determined  based  on  FE 
models  for  a  wide  range  of  finite  plate  geometries.  The  finite  width  models  included  boundary  conditions 
with  and  without  far  edge  constraints  as  an  option  to  account  for  in-plane  bending  due  to  geometric 
asymmetry.  The  resulting  solutions  are  valid  for  notch  depths  resulting  from  equivalent  hole  edge  distances 
(e/D)  from  0.0  to  49.5. 

Keywords 

Stress  Intensity  Factor  Solution 
Continuing  Damage 
Cracked  Notch 
U-Shaped  Notch 
Slotted  Notch 


1  INTRODUCTION 

When  a  crack  is  assumed  to  exist  at  a  fastener  hole  adjacent  to  the  edge  of  a  panel,  the  crack  will 
grow  to  the  edge,  resulting  in  a  slotted  edge  notch  terminated  by  the  fastener  hole.  In  most  cases, 
this  damage  is  not  sufficient  to  result  in  failure  of  the  structure.  To  continue  the  life  prediction 
process,  a  secondary  crack  is  normally  assumed  to  exist  on  the  opposite  side  of  the  hole.  This 
secondary  crack  assumption  is  the  essence  of  the  continuing  damage  concept,  as  illustrated  in 
Figure  1. 
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Figure  1 :  Typical  continuing  damage  scenario  for  a  cracked  hole  adjacent  to  the  edge  of  a  plate. 
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The  U.S.  Air  Force  defines  continuing  damage  [1]  as  the  assumption  of  a  primary  crack  (typically 
0.05  in.)  growing  from  one  side  of  a  hole  with  a  secondary  (typically  0.005  in.)  crack  growing 
simultaneously  on  the  opposite  side  of  the  hole.  Initial  through  the  thickness  cracks  are  assumed 
for  thin  skin  structures,  and  corner  cracks  are  normally  used  for  thick  sections.  Stress  intensity  (K) 
solutions  are  available  for  non-symmetric  comer  and  through  cracks  at  a  fastener  hole,  but  these 
solutions  are  only  valid  until  the  primary  crack  grows  to  the  edge  of  a  panel.  The  secondary  crack 
will  continue  to  grow,  but  the  K-solution  must  reflect  the  change  in  geometry  once  the  ligament 
between  the  hole  and  the  free  edge  is  severed. 

Closed-form  K-solutions  are  available  for  a  through  cracked  elliptical  notch  in  a  semi-infinite 
plate.  However,  it  was  determined  (based  on  notch  stress  concentration  factors)  that  these 
solutions  are  approximately  10-15%  unconservative  when  compared  to  a  U-shaped,  or  slotted 
notch,  which  is  more  representative  of  the  continuing  damage  geometry  as  shown  in  Figure  2. 


Equivalent  Notch  KT 


a  a  ld 


— U-Shapt'd  No5ch 


Figure  2:  Elliptical  vs.  u-shaped  notch  geometry  and  associated  stress  concentration  factors. 

This  paper  will  document  the  development  of  general  closed-form  K-solutions  for  corner  and 
through-the-thickness  cracks  growing  from  edge  notches  to  improve  the  accuracy  of  continuing 
damage  crack  growth  life  predictions. 


2  APPROACH 

Three  dimensional  FE  models  were  used  to  calculate  stress  intensity  factors  for  the  through  and 
corner  cracked  notch  geometries.  The  FE  analysis  program,  Abaqus  [2],  was  used  to  create  and 
analyze  all  of  the  structural  models  required  for  this  effort.  The  stress  intensity  solutions  were 
obtained  using  FRANC  3D/NG  [3]  to  insert  a  crack  in  a  local  model  surrounding  the  notch  as 
indicated  in  Figure  3. 
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Figure  3:  Illustration  of  the  global/local  FEM  approach. 

Baseline  solutions  for  a  semi-infinite  plate  were  modeled  with  a  wide  plate  as  indicated  below.  No 
units  are  shown  since  all  dimensions  are  relative. 

Baseline  model  dimensions: 

Width  =  100  *  (e  +  D/2) 

Thickness  =1.0 

Total  Plate  Height  =  500  *  (e  +  D/2) 

Crack  lengths  were  limited  to  values  that  would  result  in  (e  +  D/2  +  C)  <  10%  of  the  plate  width. 
Loading: 

Uniform  normal  pressure  along  the  top  edge 

Boundary  Conditions: 

Bottom  edge  fixed  in  the  y-direction 
Bottom,  mid-plane  node  fixed  in  the  z-direction. 

Right  edge  fixed  in  the  x-direction  (limit  in-plane  bending) 


The  finite  plate  width  effect  was  determined  by  dividing  the  FEM  results  for  various  plate  widths 
by  the  baseline  FEM  solutions  with  the  equivalent  notch  depth  and  crack  length(s). 

Solutions  for  un-constrained  in-plane  bending  were  determined  using  finite  width  FEM  solutions 
without  the  constraint  on  the  right  edge  of  the  plate. 

The  K-solution  for  all  cases  was  defined  as  a  function  of  the  remote  axial  reference  stress  Gref, 
crack  length  (x),  and  the  non-dimensional  geometry  correction  factor,  beta  (P): 


K(x)  ~  ®ref  X  /?(x) 


(1) 
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The  reference  stress  for  all  solutions  in  this  paper  is  the  gross  remote  tension  stress  applied  along 
the  top  edge  of  the  plate.  For  clarity  and  consistency  with  standard  practice,  the  crack  lengths  in 
the  plate  width  direction  were  defined  as  the  c-dimension,  and  cracks  in  the  plate  thickness 
direction  were  defined  as  the  a-dimension. 


2.1  Through-the-thickness  crack  approach 

The  2-D  (Plane  Stress)  solution  was  determined  from  the  3-D  FEM  results  using  the  standard 
relationship  between  the  plane  stress  and  plane  strain  K-solution. 

K plane  Stress  _  ■ y  ^2  ^2) 

K plane  Strain 

This  approach  was  verified  for  an  edge  crack  in  a  semi-infinite  plate  using  the  same  FEM 
dimensions,  loading,  and  boundary  conditions  used  for  the  continuing  damage  models  substituting 
the  edge  crack  length  in  place  of  the  notch  depth.  The  results  were  within  1%  of  the  well  known 
(2-D)  edge  crack  solution  [4]. 


P  =  1.122  (3) 

The  notched  solution  for  the  continuing  damage  case  is  known  when  the  crack  length  approaches 
both  zero  and  infinity.  As  the  crack  length  approaches  zero,  the  beta  value  converges  to  the  stress 
concentration  factor  for  the  u-shaped  notch  [5]  multiplied  by  the  edge  crack  correction. 


p  = 


1.122  Kt 


notch 


(4) 


As  the  crack  length  goes  to  infinity,  the  solution  converges  to  the  edge  crack  solution  in  a  semi- 
infinite  plate.  Since  the  through-the-thickness  crack  length  (c)  for  the  continuing  damage  case  is 
measured  from  the  edge  of  the  hole,  a  correction  was  applied  to  the  beta  solution  to  obtain  the 
equivalent  edge  crack  solution. 


^  1.122 

^  “  Vc/(c+W)  (5) 

The  notch  depth  (N)  was  simply  defined  as  the  hole  edge  distance  (e)  plus  the  hole  radius  (R  = 
D/2).  Wide  plate  FEM  K-solutions  were  calculated  for  several  ratios  of  hole  radius  to  notch  depth 
(R/N). 


R/N  =  0.1,  0.2,  0.333,  0.5,  and  1.0 

Crack  lengths  were  normalized  by  dividing  each  crack  length  by  the  sum  of  the  crack  length  and 
the  notch  depth.  This  allowed  the  results  for  all  notch  depths  to  be  plotted  together  so  that  data  for 
other  notch  depths  could  be  determined  by  interpolation.  It  should  be  noted  that  the  P-solution  for 
R/N=0.0  is  the  equivalent  semi-infinite  plate  edge  crack  solution. 

The  combinations  of  plate  width,  hole  diameter,  and  edge  distance  used  to  determine  the  finite 
width  correction  are  indicated  in  Table  1  and  Table  2. 
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Table  1 :  Finite  width  Geometries  for  D=0.5. 


Table  2:  Finite  width  geometries  for  D=0.25 


Hole  Dia.  =  0.25 

Plate  Width 

Offset  fe} 

0.5 

1 

1.5 

2 

4 

8 

0 

0.5 

0.75 

1 

1.5 

2 

3 

3.5 

2.2  Comer  Crack  Approach 

The  comer  cracked  solution  was  based  on  the  framework  of  the  Newman  and  Raju  closed-form  Re¬ 
solution  for  an  elliptic,  corner  cracked  hole  under  remote  tension  loading  [6]. 


K  &ref  VW  Q  Fch 


(6) 


Fch  = 


Mi + m2  q2 + m3  gy 


G1G2G3G4F(f)Fw 


(7) 


As  noted  in  reference  [6],  the  functions,  GiG2...Gn,  are  used  to  fine-tune  the  solution.  The 
functions,  Gi,  G2,  and  G3,  were  modified  as  required  to  obtain  a  close  fit  (<  5%)  to  the  FEM 
results.  Wide  plate  FEM  K-solutions  for  the  comer  cracked  hole  case  were  calculated  for  a  number 
of  R/N  values. 


R/N  =  0.1,  0.2,  0.333,  and  1.0 


The  function,  G2,  generally  accounts  for  the  unflawed  stress  distribution  in  both  crack  directions. 
G3  is  used  to  make  adjustments  for  crack  shape  (a/c)  and  percent  of  plate  thickness  penetrated 
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(a/t).  A  minor  adjustment  to  Gla  was  required  to  force  the  solution  to  converge  in  both  directions 
as  the  crack  length  goes  to  zero. 

For  the  c-dimension,  the  wide  plate  beta  values  obtained  for  the  through  crack  were  used  as  the 
basis  for  G2c. 


^2  c  - 


Pwide  Thru  Crack 


1.122 


(B) 


This  provided  a  means  to  capture  the  effect  of  a  given  notch  geometry  on  the  local  stress 
distribution  in  the  c-direction.  A  relationship  for  G2a  was  needed  that  would  allow  G3a  (required  to 
match  the  FEM  results)  to  be  within  5%  for  all  notch  depths.  The  closed-form  equation  for  G2a 
was  determined  by  iteration  using  the  FEM  results  for  each  notch  depth. 

The  solution  for  G3  was  determined  based  on  wide  plate  FEM  results  for  several  combinations  of 
crack  lengths  in  both  directions  as  indicated  in  Tables  3  through  6. 


Table  3 :  Crack  lengths  modeled  for  R/N=0. 1 . 


Table  4:  Crack  lengths  modeled  for  R/N=0.2. 


Table  5:  Crack  lengths  modeled  for  R/N=0.333. 
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Table  6:  Crack  lengths  modeled  for  R/N=1.0. 


The  range  of  crack  shapes  (a/c)  modeled  was  intended  to  cover  the  range  of  expected  values  for 
practical  structures. 

The  finite  width  correction  for  the  comer  crack  solution  was  determined  for  each  crack  dimension 
using  the  same  combination  of  plate  width  and  hole  diameter  values  that  were  used  for  the  2-D 
solution.  The  crack  shapes  for  each  model  covered  the  range  of  shape  values  used  for  the  wide 
plate  solution. 

The  general  K-solution  for  the  comer  cracked  hole  was  developed  using  FEMs  with  a  unit 
thickness  and  a  hole  radius  of  0.25.  The  effect  of  hole  radius  to  plate  thickness  (R/t)  was 
determined  based  on  additional  wide  plate  FEM  K-solutions  for  several  R/t  values. 

R/t  =  0.25,  0.50,  1.0,  2.0 

The  R/t  values  for  a  majority  of  practical  structures  fit  within  this  range.  An  additional  term,  G5, 
was  added  to  the  closed-form  solution  to  account  for  this  effect. 

3  RESULTS 

The  semi-infinite  plate  thm-crack  beta  solution  is  shown  graphically  in  Figure  4.  The  complete 
tabular  solution  is  given  in  Appendix  1 .  This  solution  serves  as  the  basis  for  the  general  through- 
the-thickness  and  comer  crack  continuing  damage  solutions.  The  tabular  solution  was  expanded 
by  interpolation  for  both  crack  length  and  notch  depth  using  FEM  results  and  the  known  limit 
solutions. 
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Figure  4:  Thru-crack  (2-D)  semi-infinite  plate  FEM  and  curve  fit  beta  solutions. 

3.1  Through-the-thickness  Crack  Results 

The  general  K-solution  for  the  thru  crack  is  expressed  as  a  function  of  the  2-D  wide  plate  beta 
solution,  the  finite  width,  and  the  in-plane  bending  constraint  factors. 

K  —  CFref  y/n  C  Pwide  Thru  Crack  (9) 

The  finite  width  correction  data  and  resulting  curve  fit  is  shown  in  Figure  5. 


Figure  5:  Thru-crack  finite  width  effect  data  and  curve  fit  results. 
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(10) 


Since  the  baseline  FEM  solutions  were  constrained  along  the  right  edge,  the  bending  correction  is 
required  to  obtain  the  K-solution  when  this  constraint  is  removed. 


X  =  (e+R+C)/W 

Figure  6:  Thru-crack  bending  correction  data  and  curve  fit  results. 


Fb  =  1  +  0.26  A  +  4.4  A2  +  11.3  A6  +  8.2  A12  +  18.5  A22  +  69.0  A39  (1 1) 


3.2  Comer  Crack  Results 

As  stated  previously,  the  corner  crack  solution  was  based  on  the  Newman  and  Raju  corner  cracked 
hole  solution  framework. 


K  =  Oref  yjna/Q  Fch 
Fch  =  k  +  M2  g)2  +  M3  g)4|  GtG2G3G4G5F^FwFh 


(12) 

(13) 


The  majority  of  the  parameters  were  used  without  changes  and  are  given  here  as  a  convenience  for 
the  reader. 


Fora/c  <  1 


M1  =  1.13  -  0.09  (a/c) 

(14) 

0-89  _  . 

Mo  =  7 - —  -  0.54 

z  (0.2+a/c) 

(15) 

9 

Approved  for  public  release;  distribution  unlimited. 


M3  =  °-S-<S^+14(1-«A')24  (16) 

Glc  =  1  +  (0.1  +  0.35  (a/t)2)  (17) 

F<j>c  =  yja/c  (18) 

F<pa  =  1.0  (19) 

Q  =  1.0  +  1.464  (a/c)1-65  (20) 

Fora/c  >  1 

M1  =  JcJa  (1.0  +  0.04  (c/a))  (21) 

M2  =  0.2  (c/a)4  (22) 

M3  =  -  0.11  (c/a)4  (23) 

Clc  =  1  +  (0.1  +  0.35  (c/a)  (a/t)2)  (24) 

F0C  =  1.0  (25) 

F(pa  =  4cja  (26) 

Q  =  1.0  +  1.464  (c/a)165  (27) 


The  remaining  parameters  were  modified  as  required  to  provide  the  best  fit  to  the  FEM  results. 
Gla  =  1.0  + 0.0699  e-150^ 


^2  C  - 


Pwide  Thru  Crack 


1.122 


G2 a  ~  G 


2c 


0.985  + 


1.31-0.985 


—  0.323— 7-fjrpr 
1+e  V  (C+N), 


1/0.2 


(28) 

(29) 

(30) 
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1.2 


0  0.2  0.4  0.6  0.8  1 


♦  A/C=1.0 
■  A/C=0.75 
▲  A/C=0.5 
X  A/C=2.0 
X  A/C=3.0 

- Fit,  A/C=1.0 

- Fit,  A/C=0.75 

- Fit,  A/C=0.5 

- Fit,  A/C=2.0 

- Fit,  A/C=3.0 


a/t 


Figure  7:  Comer  crack  G3c  FEM  data  and  curve  fit  results. 


Fora/c  <  1 

G3c  =  (0.985  -  0.087  (a/c)425)  +  0.29  (a/t)  -  0.11(a/t)425 
0.05  (a/t)6 


Fora/c  >  1 


G3c  =  (1.5585  -  0.66  (a/c)01)  +  0.29  (a/t)  -  0.11(a/t)425 
0.05  (a/t)6 


0  0.2  0.4  0.6  0.8  1 


♦  A/C=1.0 
■  A/C=0.75 
▲  A/C=05 
X  A/C=2.0 
X  A/C=3.0 

- Fit,  A/C=1.0 

- Fit,  A/C=0.75 

- Fit,  A/C=0.5 

- Fit,  A/C=2.0 

- Fit,  A/C=3.0 


a/t 


Figure  8:  Comer  Crack  G3a  FEM  data  and  curve  fit  results. 
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Fora/c  <  1 


63 a  ~  0-735  + 


t 


[l.625-0.4(a/c)225](a/t) 


175+(a/t)(1-32-°-22(“/c): 


y 


+  0.625  (a/t)e 


(33) 


Fora/c  >  1 


63 a  ~  0.735  + 


/[i.24-0.0i5(a/c)ls](a/t)\  ft„cf  8 
\0.175+(a/t)(112_002^a/c))/  +  °‘625  (a/t) 


(34) 


•  R/T  =0.25 

♦  R/T  =  0.5 
▲  R/T  =1.0 
■  R/T  =2.0 

- Fit,  R/T  =  0.25 

- Fit,  R/T  =  0.5 

- Fit,  R/T  =  1.0 

- Fit,  R/T  =  2.0 


Figure  9:  Comer  Crack  G5c  FEM  data  and  curve  fit  results, 
Gse  =  1  +  (  1S4W;);‘  -  0.035)  (1  -  {a/0f 

bc  Vl+10.445  (R  t)2A  )  v  v  J  J 


25 


(35) 


0.2  0.4  0.6 

a/t 


•  R/T  =0.25 

♦  R/T  =  0.5 
▲  R/T  =  1.0 
■  R/T  =2.0 

- Fit,  R/T  =  0.25 

- Fit,  R/T  =  0.5 

- Fit,  R/T  =  1.0 

- Fit,  R/T  =  2.0 


Figure  10:  Comer  Crack  G5a  FEM  data  and  curve  fit  results. 
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c~  =  °'5 + (JT, :.«/.>))  +  ((“  w,a)  (W«  -  °-25)  <»/«2) 


(36) 


♦  A/C=3_0 
■  A/C=2_0 
A  A/C=1.0 

*  A/C-0.5 
- Fit 


0  0.2  0.4  0.6  0.8  1 

X  =  (e+R+C)/W 

Figure  1 1 :  Comer  Crack  Fwc  data  and  curve  fit  results. 


F  = 

1  wc 


(37) 


X  =  (e+R+C)/W 


Figure  12:  Comer  Crack  Fwa  data  and  curve  fit  results. 


F  = 

1  wa 


(38) 
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The  in-plane  bending  correction  is  based  on  the  curve-fit  for  the  through  crack  case  (Fb,  Section 
3.1)  where  the  function,  X,  is  defined  for  each  crack  direction. 


Fbc  ~  Fb(Aa)  +  (a/t)2  (FbUc)  -  Fb(Xa)) 

(39) 

Fba  =  Fb(A la)  +  (Fbac )  -  FbUa )) 

(40) 

Ac  =  (e  +  R  +  c)  /  W 

(41) 

Aa  =  (e  +  R)/W 

(42) 

♦  C-Dimension 

•  A-Dimension 


C/(W-Notch  Depth) 


Figure  13:  Corner  Crack  in-plane  bending  correction  error  vs.  c-length  to  un-cracked  net  section 
ratio. 


4  VERIFICATION 

The  2-D  tabular  semi-infinite  plate  beta  solution  matrix  was  expanded  to  include  a  wide  range  of 
notch  depths  and  a  large  number  of  crack  lengths.  The  values  were  selected  so  that  reasonably 
accurate  (within  5%)  solution  for  a  given  geometry  could  be  obtained  using  linear  interpolation. 

Since  the  proposed  through-the-thickness  and  comer  cracked  solutions  were  developed  by  curve 
fitting  the  baseline  FEM  results,  good  agreement  for  all  baseline  notch  depths  would  be  expected. 
The  solutions  were  verified  for  a  notch  depth  that  was  significantly  deeper  (R/N=0.03)  than  all  of 
the  baseline  FEM  models.  A  limited  number  of  crack  lengths  were  analyzed  since  CPU  mn  times 
were  approximately  12  hours  for  each  case  due  to  the  size  of  the  models. 
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C/(C+M) 

Figure  14:  Wide  plate  through-the-thickness  crack  beta  solution  verification  for  R/N  =  0.03. 


Table  7:  Comer  crack  FEM  geometry  used  for  K-solution  verification 


Plate  Width 

Hole  Dia. 

Thickness 

R/N 

50.0 

0.5 

1.0 

0.030 

Table  8:  Comer  crack  verification  results. 


c 

a 

a/c 

Kc  (Fit) 

Kc  (FEM) 

%  Diff 

Ka  (Fit) 

Ka  (FEM) 

%  Diff 

0.025 

0.025 

1.000 

2.202 

2.226 

-1.167 

2.466 

2.500 

-0.560 

0.050 

0.050 

1.000 

2.750 

2.767 

-1.326 

3.516 

3.420 

2.607 

0.100 

0.100 

1.000 

3.133 

3.271 

-2.365 

4.756 

4.574 

4.023 

0.100 

0.300 

3.000 

4.316 

4.536 

-4.632 

4.656 

4.443 

4.734 

0.250 

0.250 

1.000 

3.774 

3.603 

-0.763 

6,561 

6,342 

3.453 

0.250 

0.500 

2.000 

4.347 

4.360 

-0.663 

6.417 

6,237 

1.306 

0.500 

0.500 

1.000 

4.476 

4.546 

-1.563 

7.643 

7.762 

-1.533 

0.500 

0.250 

0.500 

2.672 

2.661 

-0.336 

7.367 

7.353 

0.103 

0.300 

0.300 

1.000 

5.636 

5.745 

-0.653 

3.720 

10.255 

-5.217 

5  DISCUSSION  AND  CONCLUSIONS 

Baseline  FEM  data  were  limited  to  R/N  values  from  0.1  to  1.0.  This  corresponds  to  an  edge 
distance  (e/D)  range  of  4.5  to  0.0.  While  this  would  be  sufficient  for  most  practical  cases,  a  Re¬ 
solution  is  often  required  for  higher  e/Ds  when  a  secondary  crack  grows  along  a  row  of  holes  prior 
to  final  fracture.  Extrapolating  the  2-D  semi-infinite  plate  results  to  an  R/N  of  0.01  allows  a 
solution  to  be  calculated  for  edge  distances  as  high  as  49.5.  Although  FEMs  for  R/N=0.01  were 
not  attempted  as  part  of  this  effort,  there  was  very  good  agreement  between  the  FEM  results  and 
the  extrapolated  beta-values  for  R/N=0.03.  In  addition,  the  extrapolation  was  well  bounded  by  the 
known  equivalent  solution  for  the  edge  crack  (R/N=0.0).  The  equivalent  continuing  damage  beta 
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solution  for  an  edge  crack  must  be  corrected  since  the  crack  length  is  measured  from  the  edge  of 
the  notch.  The  2-D  semi-infinite  plate  beta  table  indicates  that  the  solution  is  nearly  converged  to 
the  equivalent  edge  crack  solution  for  notch  depths  beyond  R/N=0.01. 

The  in-plane  bending  correction  for  the  comer  crack  solution  is  complicated  by  the  out-of-plane 
bending  effect  of  the  corner  crack  geometry.  The  proposed  corner  crack  bending  correction  is  in 
good  agreement  (<  6%)  with  the  FEM  results  when  the  crack  length  in  the  c-direction  is  less  than 
20%  of  the  un-cracked  net  section  width.  There  are  few  practical  cases  that  would  be  affected  by 
this  limitation. 

The  proposed  closed-form  solutions  are  offered  as  a  reasonably  accurate  alternative  to  the  use  of 
detailed  FEM  K-solutions  for  continuing  damage  crack  growth  life  prediction. 
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Appendix  1:  Thru-Crack  (2-D)  Semi-Infinite  Plate  Beta  Table 


2-D  Semi-Infinite  Plate  Betas 


C/(C+N) 

K/N 

0 

0.01 

0-02 

0,04 

0.0625 

04 

045 

0.2 

0,25 

033333 

0,5 

0-75 

1 

0.000 

N/A 

25.9205 

18.6830 

13.5774 

11,1368 

8.9768 

7.5068 

6,6428 

6.0425 

5,3580 

4.5389 

3.8600 

3.4392 

0.003 

20.4848 

18.1250 

14.8275 

11.9950 

10.2450 

8.5050 

7.2365 

6.4478 

5.8677 

5,2290 

4.4753 

3.8185 

3.4108 

0.00625 

14.1923 

13  9750 

12.4500 

10.7300 

9.4750 

8  0625 

6  9678 

6.2525 

5.6950 

5.1050 

4.4127 

3  7750 

3.3820 

0.009 

11.8269 

11.8000 

10,9100 

9,3650 

8.8948 

7.7175 

6,7545 

6,0965 

5.5635 

5.0095 

4,3609 

3.7393 

3.3583 

0.0125 

10.0355 

10.0355 

9.5000 

9.0000 

S.2500 

7.3220 

6.5000 

5,9100 

5.4100 

4.9000 

4.2965 

3,6950 

3,3310 

Q.Q1875 

8.1939 

8.1939 

8.0350 

7.7150 

7.3100 

6.6985 

6.0750 

5,5900 

5.1500 

4.7225 

4.1855 

3.6180 

3.2822 

0.0250 

7.0962 

7,0962 

7.0962 

6,3500 

6.6000 

6.1920 

5.7000 

5,3050 

4,9250 

4,5750 

4,0800 

3.5450 

3.2375 

0.0375 

5.7940 

5.7940 

5.7940 

5.7600 

5.6700 

5.4500 

5.1450 

4.8650 

4.6000 

4,3000 

3.9000 

3.4200 

3.1535 

0.050 

5.0177 

5.0177 

5.0177 

5.0177 

49500 

4.8400 

4.6500 

4.4800 

4.3125 

4.0700 

3.7255 

3,3150 

3,0738 

0.0625 

4,4880 

4.4880 

4.4880 

4.4880 

4.4880 

4.3750 

4.2800 

4.1700 

4.060C 

3.8595 

3.5750 

3  2200 

2  9980 

0.075 

4,0970 

4,0970 

4.0970 

4,0970 

4,0970 

4.0350 

3.9820 

3,9220 

3.8350 

3,6750 

3,4320 

3.1300 

2.9315 

0.0875 

3.7931 

3.7931 

3.7931 

3.7931 

3.7931 

3.7755 

3.7415 

3.7000 

3.6375 

3,5100 

3.3050 

3.0450 

2.8676 

0.100 

3.5481 

3,5481 

3.5481 

3,5481 

3.5481 

3.5481 

3.5325 

3.5050 

3.4600 

3.3575 

3.1880 

2.9650 

2.8040 

0.125 

3.1735 

3.1735 

3.1735 

3,1735 

3.1735 

3.173S 

3.1680 

3.1625 

3.1300 

3,0750 

2.9650 

2.8145 

2.6800 

0.150 

2.8970 

2,8970 

2.8970 

2,8970 

2,8970 

2.8970 

2,8970 

2,8950 

2.6900 

2,8600 

2.7850 

2.6675 

2.5600 

0.175 

2.6821 

2.6821 

2.6821 

2.6821 

2.6821 

2.6821 

2.6821 

2.6821 

2.6800 

2.6728 

2.6215 

2  5400 

2.4495 

0.200 

2.5089 

2.5089 

2.5089 

2,5089 

2.5089 

2.5089 

2.5089 

2,5089 

2.5089 

2.5070 

2,4750 

2.4150 

2.3490 

0.250 

2.2440 

2.2440 

2.2440 

2.2440 

2,2440 

2,2440 

2  2440 

2.2440 

2.2440 

2.2440 

2.2375 

2.2100 

2.1739 

0,300 

2.0485 

2,0485 

2.0485 

2,0485 

2.0485 

2.0485 

2.0485 

2.0485 

2.0485 

2,0485 

2.0420 

2.0300 

2.0185 

0.350 

1.8965 

1,8965 

1.8965 

1,3965 

1.8965 

1.8965 

1.8965 

1,8965 

1.6965 

1,8965 

1,8965 

1.8920 

1.8870 

0,400 

1.7740 

1.7740 

1.7740 

1.7740 

1.7740 

1.7740 

1.7740 

1.7740 

1.7740 

1.7740 

1.7740 

1.7738 

1.7730 

0.450 

1,6726 

1.6726 

1.6726 

1.6726 

1.6726 

1.6726 

1,6726 

1,6726 

1.6726 

1.6726 

1,6726 

1.6726 

1.6695 

0.500 

1,5867 

1.5867 

1.5867 

1.5867 

1  5867 

1  5867 

1.5867 

1.5867 

1.5867 

1.5867 

1.5867 

1.5867 

1.5867 

0.550 

1,5129 

1.5129 

1.5129 

1,5129 

1.5129 

1.5129 

1,5129 

1,5129 

1.5129 

1.5129 

1,5129 

1.5129 

1.5129 

0.600 

1.4485 

1.4485 

1.4485 

1.4485 

1.4485 

1.4485 

1.4485 

1.4485 

1.4485 

1,4485 

1.4485 

1.4485 

1.4485 

0.700 

1.3410 

1.3410 

1.3410 

1.3410 

1.3410 

1  3410 

1  3410 

1.3410 

1.3410 

1.3410 

1.3410 

1.3410 

1  3410 

0.800 

1,2544 

1,2544 

1.2544 

1,2544 

1,2544 

1.2544 

1,2544 

1,2544 

1.2544 

1,2544 

1,2544 

1.2544 

1.2544 

0,900 

1.1827 

1.1827 

1.1827 

1,1827 

1.1827 

1.1827 

1.1827 

1.1827 

1.1827 

1,1827 

1.1827 

1.1827 

1.1827 

1.000 

1.1220 

1.1220 

1.1220 

1.1220 

1.1220 

1.1220 

1.1220 

1.1220 

1.1220 

1.1220 

1.1220 

1.1220 

1.1220 

The  shaded  area  indicates  that 
solution  as  noted  in  Section  2.1. 


the  beta  solution  has  converged  to  the  equivalent  edge  crack 
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